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ABSTRACT 

We have mapped 20 molecular lines in the Central Molecular Zone (CMZ) around the 
Galactic Centre, emitting from 85.3 to 93.3 GHz. This work used the 22-m Mopra radio 
telescope in Australia, equipped with the 8-GHz bandwidth UNSW-MOPS digital filter bank, 
obtaining ^ 2 km s^^ spectral and ^ 40 arcsec spatial resolution. The lines measured include 
emission from the C-C3H2, CH3CCH, HOCO+, SO, H^^CN, Hi3CO+, SO, H^^NC, C2H, 
HNCO, HCN, HCO+, HNC, HC3N, "CS and N2H+ molecules. The area covered is Galactic 
longitude —0.7 to 1.8 deg. and latitude —0.3 to 0.2 deg., including the bright dust cores around 
Sgr A, Sgr B2, Sgr C and Gl. 6-0. 025. We present images from this study and conduct a 
principal component analysis on the integrated emission from the brightest 8 lines. This is 
dominated by the first component, showing that the large-scale distribution of all molecules 
are very similar. We examine the line ratios and optical depths in selected apertures around 
the bright dust cores, as well as for the complete mapped region of the CMZ. We highlight the 
behaviour of the bright HCN, HNC and HCO+ line emission, together with that from the ^'^C 
isotopologues of these species, and compare the behaviour with that found in extra-galactic 
sources where the emission is unresolved spatially. We also find that the isotopologue line 
ratios (e.g. HCO"'"/H^'^CO"'") rise significantly with increasing red-shifted velocity in some 
locations. Line luminosities are also calculated and compared to that of CO, as well as to line 
luminosities determined for external galaxies. 

Key words: ISM:molecules - radio lines:ISM - ISM:kinematics and dynamics. 



1 INTRODUCTION 

The Central Molecular Zone (CMZ) is th e region within about 300 
parsec (2 degrees) of the Galactic Centre I IMorris &Serabvnll 19961) 
with a strong concentration of molecular g as. The molecular distri- 
bution is most clearly traced by CO (e.g. Ballv et al. 'I987','l988'; 
Bitran etal. 1997: O ka et al. 1998; Dahmen et al. 1997; Oka et al. 
20071 ; iMartin et al ll2004) . which shows high gas densities, large 
shear velocities in the Galactic Centre potential well with non- 
circular motions, large internal velocity dispersions within the 
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clouds and high gas temperatures. Th e molecular mass in th e CMZ 
is estimated as (3^i) x 10^ Mq bv lOahmen et al.l ( 1 19981) from a 
range of tracers. 

Some of the most prominent features jPalmer & Goss|[T996h 
of the CMZ are Sagittarius A (Sgr A) the Galactic Centre cloud, 
Sagittarius B2 (Sgr B2) and Sagittarius C (Sgr C) , named from 
the radio and near infrared peaks (P iddington & MinnetJ 1 195 ll ; 
lLeaueuxlll962l ; lHoffmatm. Frederick & Emery.. 197 ll) . 

The CMZ is asymmetric about the Galactic Centre, with the 
molecular distribution biased towards positive longitudes, with 
Sgr B2 a major feature at / = +0.68° . The three-dim ensional shape 
has been modelled as an elongated structure e.g. bv lSawada et al.l 
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( I2004I) using constraints from the CO emission and OH absorption, 
witii tiie positive longitude end closer to us. 

The CMZ is prominent in far-infrared and sub-millimetre 
dust emission, e.g. JCM T SCUBA observations at 850 and 450 
um jPierce-Price et alj I2OOO). APEX ATLASGAL at 870 /xm 
dSchuUer et alj|2009|) , Bolocam at 1100 and 350 ^.m ( iBallv et al.l 
I2OIOI) and HiGAlFI at 70, 170, 250, 350 and 500 ^m. This 
traces the T ~20 K du st associated with the molecular gas with 
IPierce-Price et al.l 1 I2OOOI) estimating the total gas mass in the CMZ 

as (5.3 ± 1.0) X lO'^ Mq from the dust emission^ 

The CMZ is a site of recent star formation jYusef-Zadeh et all 

I2OO9I) and very high stellar densities. Because of the very high ex- 
tinction and crowding, the stellar objects are best studied in the 
infrared and wit h high resolution, e.g . see Spitzer studies at 24 ^m 
jHinz et alj|2009i) and 3.6 to 8.0 /^m jRamirez et a"i]|2008l) . 

The molecular medium of the CMZ has a surpr isingly rich 
chemistry. The Sgr B2 region is well-studied (e.g. Ijones et aTl 
|2008a, 20nj, and references therein) and known to have rich or- 
anic chemistry, particularly the Sgr B2 (N) and (M) hot cores 
Herbst & van DishoecM 120091) . However, complex molecules are 
extended over the larger CMZ area as s hown by widesprea d 
CH3OH (methanol) at 834 MHz detected bv lOottlieb etaL ( 1979), 
HNCO (isocyanic acid) at 110 GHz imaged by (Dahmen et al. 
j 19971) a nd complex organic molecu les detected in individual 
clouds bv lReguena-Torres et al. Il l2008h . The chemistry in the CMZ 
is similar to that in hot molecular cores (HMCs, hot star-forming 
dense cores withi n giant molecular cloud com plexes) but over a 
much larger scale jRequena-Torres et al.ll20od) . 

Other molecule s ima ged in the CMZ include SiO 
jMar tin-Pintado et al.' 'l997 ), wh ich tra ces extensive shocks. 



CS iTsuboi, Handa & Ukit j|l999l) . NH3 jNaga vama et al. 2 00SD, 
HCN jjackson et al. 199 C^HCO+ and H^^CO+ dRiguelme eTS 
l2010ah and OH tBovce & Coheij|l994l) in absorption against the 

continuum. 

It was argued bv^Htittem eisteJ ( fl998l) that "large scale surveys 
in as many molecules, isotopomers and transitions as possible are 
essential to understand the structure of the CMZ" preferably with 
resolution < 2 pc (< 50 arcsec). The reason is that the distribu- 
tion of a particular line depends on the excitation, optical depth and 
abundance. The '^^CO 1-0 line, in particular, does not trace the 
density structure well. To obtain the density and temperature distri- 
bution of the gas, multilevel studies, isotope studies and surveys in 
rarer, weaker molecular species are required. 

We aim, in part, to fill this need. We present here observa- 
tions of the CMZ of 20 lines in the 3-mm band (including two 
hydrogen recombination lines), with resolution around 40 arcsec. 
Some preliminary results from this pro ject have been presented in 
Ijones, Burton. & CunninghamI ( l2008bb and ljones. Burton. & Lowel 
( l2008d) . 

This is part of a larger project in which we are also observ- 
ing the CMZ area in multiple lines in the 7-mm band, with resolu- 
tion around 65 arcsec. Additionally we have imaged the CMZ area 
in multiple lines in the 12-mm band, with resolution around 120 
arcsec, as part of the wide r H2O southern Galactic Plane Survey 
(HOPS) ( Walsh et al.ll2008ll201 ih . 

Due to the particular chemical complexity, and line strength, 
of the Sgr B2 region, we have multi-line spectral images of this 
smaller sub-area of the CMZ over the full 3-mm band from 82 



to 114 GHz, and 7-m m band from 30 t o 50 GHz, presented in 
Ijones et al] l l2008ah and ljones et alj ( l201lh respectively. 

We use the distance to the Galactic Centre of 8.0 kpc 
( iReid et alj2009h . so 1 arcmin corresponds to 2.3 pc. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were made with the 22-m Mopra radio tele- 
scope, of the Australia Telescope National Facility (ATNF) using 
the MOPS digital filterbank. The Mopra MMIC (Monolithic Mi- 
crowave Integrated Circuit) receiver has a bandwidth of 8 GHz, 
and the MOPS backend can cover the full 8-GHz range simultane- 
ously in the broad band mode. This gives four 2.2-GHz sub-bands 
each with 8192 channels of 0.27 MHz. The lines in the CMZ are 
broad, so that the 0.27 MHz channels, corresponding to around 
0.9 km s^^, are quite adequate. 

The Mopra receiver covers the range 77 to 1 17 GHz in the 3- 
mm band. We chose the tuning centred at 89.3 GHz, to cover the 
range 85.3 to 93.3 GHz, giving the spectral lines summarised in 
Table[T] This range was chosen to include the strong lines of HCN, 
HCO^ and HNC, which are some of the strongest lines in the 3-mm 
band (after ^^CO and ^'^CO), plus a good range of other lines. 

The area was observed with on-the-fly (OTP) mapping, in 
blocks of 5 . 1 X 5. 1 arcmin^, in a similar way to that described in 
Ijones et al] ( l2008ah . We used position switching for bandpass cali- 
bration with the off-source reference position observed before each 
5 arcmin long source scan. The reference position ((a, (5)j2ooo = 
17''51™03.^ 6, -28°22'47", or I = 1.093 deg., b = -0.735 deg.), 
was carefully selected. It is hard to balance the compromise of po- 
sitions which are relatively emission free, in the Galactic Centre 
area, without too large a position offset for good spectral baselines. 
The OTP maps used scan rate 4 arcsec/second, with 12 arcsec line 
spacing, taking a little over an hour per block. We made pointing 
observations of SiO maser positions (AH Sco or VX Sgr), before 
every map, to correct the pointing to within about 10 arcsec ac- 
curacy. The system temperature was calibrated with a continuous 
noise diode, and ambient-temperature load (paddle) about every 30 
minutes. The system temperatures were mean 210 to 225 K across 
the 85 to 93 GHz band, with standard deviation 50 K, depending 
on elevation and weather conditions, and extremes 155 to 410 K. 

The 5x5 arcmin^ blocks were observed twice each with 
Galactic latitude and longitude scan directions, slightly offset, to re- 
duce scanning direction stripes, and to improve the signal to noise. 
The overall 2.5 x 0.5 deg'^ area was covered with a 30 x 6 grid of 
blocks, separated by 5 arcmin steps, making 360 OTP observations 
required. In practice close to 400 OTP observations were used, in- 
cluding some areas with OTP maps stopped by weather or other 
problems, and re-observed. The area covered is between -0.72 to 
1.80 degrees Galactic longitude, and -0.30 to 0.22 degrees Galactic 
latitude. 

The observations were spread over three southern winter sea- 
sons, in 2007, 2008 and 2009. Although the blocks worst affected 
by poor weather were re-observed, the range of conditions dur- 
ing the different observations means that there are inevitably some 
blocks with greater noise level than others. 

The OTP data were turned into PITS data cubes with the LIVE- 
DATA and GRIDZILLA package^ The raw spectra in RPPIT^ for- 
mat, were bandpass corrected and calibrated using the off-source 



^ https://hi-gal.ifsi-roma.inaf.it/higal/ 



^ http://www.atnf.csiro.au/people/mcalabre/livedata.html 
^ http ■J/vi'ww. atnf .csiro . au/computing/soft ware/rpfits .html 
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reference spectra with LIVEDATA, then a robust second order poly- 
nomial fitted to the b aseline an d subtracted, with the data output 
formatted as SDFITS ( [Garwood 2000) spectra. These spectra were 
then gridded into datacubes using GRIDZILLA, with a median fil- 
ter for the interpolation, and combination of data over-sampled in 
position. The median was used, as this is more robust to the out- 
liers caused by bad data. The scripts used for gridding allowed the 
lines (Table [TJ to be specfied, with their rest frequencies, so the 
GRIDZILLA output was FITS cubes with velocity coordinates, com- 
bining data for the whole mapped area. 

The FITS cubes were then read into the MIRIAD package for 
further processing and analysis. In particular, as the emission is typ- 
ically of low surface brightness, the data were smoothed in velocity, 
with a 7-point banning kernel, to make a version of the data cubes 
with improved surface brightness sensitivit3Q. This gives around 

1.07 MHz, or 3.6 km s~^ effective spectral resolution, which we 
use as a reduced size data cube by droppingevery second origi- 
nal pixel to make a Nyquist sampled versior|j with 0.54 MHz or 

1.8 km s~^ pixels. This is still quite adequate spectral resoution 
for the broad lines in the CMZ, which are > 10 km s~^ wide, but 
for the rare narrow spectral feature, we use the data cubes in the 
original 0.27 MHz or 0.9 km s^^ pixel versions. 

The resolution of the Mopra be am varies between 36 arcsec at 
86 GHz and 33 arcsec at 115 GHz iLadd et alj|2005l) . so the res- 
olution in the final data is around 39 arcsec, after the effect of the 
median filter convolution in the gridzilla interpolation. The main 
beam efficiency of Mopra varies between 0.49 at 86 GHz and 0.44 
at 100 GHz (Ladd et_ al., 2005 , ). How ever, the Mopra beam has sub- 
stantial sidelobes jLadd et alj2005l) so that the extended beam effi- 
ciency, appropriate for the extended emission we are mostly study- 
ing here, is quite different (0.65 at 86 GHz). Since we are largely 
concerned in this paper with the spatial and velocity structure, we 
have mostly left the intensities in the T^ scale, without correction 
for the beam efficiency onto the Tm b scale, except for some of the 
quantitative analysis sections. 



3 RESULTS 

3.1 Line data cubes 

We have data cubes for the twent)lf| strongest lines in the 85.3 
to 93.3 GHz range, as listed in Table [T] The line identifica- 
tions with rest frequencies and transitions are taken from the on- 
line NIST cataloguj^ of lines known in the interstellar medium 
jLovas & Dragoseil 20041) and the splatalogue compilatior[f|. 

The root-mean-square (RMS) noise level and peak brightness 
are listed in Table |2l The RMS noise is around 50 mK, except for 
the two lines of C2H at 87.32 and 87.40 GHz which are close to the 
edge of the sub-band at 87.3 GHz, where there is poor data. 

The velocity range listed in Table [2] indicates the range of sig- 
nificant emission detected in the data cubes, and was the range over 

* The sensitivity of brightness temperature in K is improved by increasing 
the channel bandwidth. 

^ The 7-point banning smoothing gives FWHM four pixels, so dropping 
alternate pixels gives two pixels per FWHM. 

^ As the observations included the full 8-GHz range, there are other weaker 
lines detected, mostly in Sgr B2, but as they are quite weak and not extended 
over the l arger area of the CMZ, we do not consider them further here. See 
Table 3 of ljones et aj] f2008a) for a listing of these Sgr B2 lines. 
http://physics.nist.gov/PhysRefData/Micro/Html/contents.html 

* http://www.splatalogue.net/ 



Table 1. The lines imaged here in the 85.3 to 93.3 GHz range. The multiple 
components listed here are, in general, blended due the large velocity widths 
in the CMZ area. The observations covered the whole 85.3 to 93.3 GHz 
range, s o do include many more weaker lines, particulai'ly in Sgr B2, as 
listed in I Jones et all )2008j) . 
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which we integrated to get images of the total line emission. We 
show the spectrum of HCN averaged over the whole CMZ area in 
Fig. [T] which has strong emission, seen over the largest velocity 
range. The range of velocities is large, due to the large line widths 
in the CMZ and the large velocity gradient across the CMZ in the 
deep potential well of the Galactic Centre. This makes the inte- 
grated emission image quite sensitive to low level baselevel offsets, 
which while small in terms of brightness (e.g. of order 10 mK) be- 
comes significant (e.g. of order K km s~^) when integrated over the 
velocity range (order 100 km s^^). See section [34] for discussion 
of the integrated emission. 
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Table 2. Statistics of the data cubes. The RMS noise values are from the hanning smoothed data with 1.07 MHz pixels: the original 0.27 MHz pixel data have 
noise a factor of 2 greater Both the RMS noise and peak brightness temperature are in TJ . The velocity range is that with line emission significantly above 
the noise level, and is the range used to integrate to get total line emission images. We also quote the position and velocity of the peak pixel. 
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* = line wing coincides with bad data at sub-band edge. 



noise is reduced in terms of brightness, so that weak emission is 
then clearer above the reduced noise. 

We also note that there are some artifacts in some of these peak 
brightness images, notably stripes in the latitude scanning direction 
due to data taken in the poorest weather, around I — 1.4 deg., b — 
0.2 deg. and / = 359.5 deg., 6 = -0.1 deg. 

The peak emission images do not show the line absorption. 
There is absorption in some stronger lines (e.g HCN, HCO^ and 
HNC) at the strong continuum of Sgr B2 jjones et al.ll2008ij) and 
over a wider area, see subsection l3.3l and Fig.|7] However, the peak 
emission is mostly not sensitive to this absorption, as the velocities 
of peak emission are offset from that of the absorbing clouds. 

The lines of HCN, HCO+ and HNC (Fig. Hi are strong and 
detected over most of the imaged area. The ^"'C isotopologues, 
H^^CN, H"C0+ and HN^^C (Fig. [3) are detected in the densest 
cores, mostly around Sgr B2 and Sgr A. The ratio between these 
and the HCN, HCO+ and HNC lines shows that the latter are opti- 
cally thick in the densest cores (subsection |3.6t . This optical depth 
effect explains why the position of the brightest peak of HCN, and 
the velocity of the brightest peak of H^'^CN are different to that of 
other lines. 

Other strong lines that are detected over a large fraction of the 
CMZ are HNCO, N2H+, HC3N (Fig.[2j, CH3CN and SiO (Fig.^. 
The HNCO and HC3N peaks are particularly strong around Sgr B2 
(Tablellll. 

Weaker lines that are detected mostly around Sgr B2 and Sgr A 
are HOCO+, ^^CS, c-C3H2,CH3CCH and C2H (Fig.H. The data 
for the two lines of the C2H molecule show noise and greater arti- 
facts due to the line frequencies being at the edge of the sub-band, 
giving poor sensitivity and baseline stability. 

The line of SO (Fig. |4j is concentrated at Sgr B2, as are 
the two radio recombination lines, H 41 a and H 42 q (Fig. [Sj. 
These recombination lines indicate significant free-free emission 

© 201 1 RAS, MNRAS O00.mi27l 



HCN mean spectrum 




-300 -200 -100 100 200 300 



Velocity / km/s 

Figure 1. Spectrum of HCN 1-0 averaged over the whole 2.5 X 0.5 deg^ 
CMZ area observed. 

3.2 Peak emission images 

We show the 2-dimensional distribution of the line emission (as 
listed in Table [T] with quantum numbers) using images of the peak 
emission in Figs|2]to|5] These images are much more robust to the 
effect of low level baseline problems than the integrated emission 
images. They also have the advantage of showing fine structure, 
being less affected by multiple velocity components. We use the 
peak brightness within the velocity ranges listed in Table |2] The 
grey-scale displays start at the 3 a brightness level, so that the dis- 
play suppresses areas without significant line emission. However, it 
should be noted that there is emission outside the areas highlighted 
in these plots, below the 3 a level for individual pixels, which is 
found to be significant when the data are integrated over larger ar- 
eas than the beam. That is, when averaging over larger areas, the 



3 -mm spectral imaging of the CMZ 5 



HCN 88,63 GHz 




NgH+ 93.17 GHz 




1.5 1 0.5 359.5 

GLON (degrees) 



Figure 2. Peak brightness images, for the lines of HCN, HCO+, HNC, HNCO, N2H+ and HC3N. These are the six brightest lines we observed. The grey-scale 
is peak brightness as TJ in K here and in Figs[3]to|5] 
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Figure 3. Peak brightness images, for the lines of H^^CN, H13cO+, HNI^C, CH3CN, SiO and HOCO+. 
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at 3 mm at Sgr B2, as the cascade through Rydberg states of H 
of recombinitig electrons occurs under similar conditions to the 
thermal bremsstrahlung of the free electrons in the ionised gas 
dWilson. Rohlfs. Huttemeiste3l2009l) . The free-free emission com- 
ponent of the Sgr B2 spectrum at 3 m m is also shown by the overall 
spectral energy distribution (SED) in ljones et alj ( l201lh . 

The area around Sgr B2 has interesting differences in the 
line distributions, which are not well shown at the scale of Figs 
2]to[5l but which are shown and discussed in the previous paper 
Jones et alj fcoOSal) . 

3.3 Velocity structure 

The line emission in the CMZ spans a large range of velocities 
(Fig. [T), which, as we note in subsection 13.11 is due to both the 
large gradient in velocity across the CMZ in the deep potential well 
of the Galactic Centre and large line widths in the CMZ. We show 
in Fig.[6]images of the HCN line emission as a function of velocity, 
integrated over 36 km s^^ ranges. This demonstrates the large ve- 
locity gradient across the Galactic Centre, with negative velocities 
at more negative longitudes, and positive velocities at more positive 
longitudes. 

The molecular line emission in the CMZ is not symmetrically 
distributed around the Galactic Centre, so the strongest emission 
appears at positive longitudes (dominated by Sgr B2 at around lon- 
gitude 0.7 degrees) giving an asymmetric distribution in velocity 
(about Vlsr = km s^^) in Fig.[6]and the profile in Fig.[T] The 
area we have observed to cover the strongest emission is longitude 
—0.7 to 1.8 deg., so there is also some observational selection in 
this bias to positive longitudes and velocities. 

Another way to show the velocity structure is with a position- 
velocity diagram, as in Fig. |7] for the peak brightness of HCN as 
a function of velocity and longitude, within three latitude ranges 
10 arcmin thick covering -0.29 to 0.21 deg. This shows both the 
velocity gradient across the region, as well as large velocity widths 
at some locations. 

Also shown in Fig. |7] are absorption features at velocities 



around —52, —28 and —3 km s~^, which are also seen in the in- 
tegrated emission spectrum (FiglTJ. These absorption features are 
clear in the channel images of the HCN data cube, but are not ob- 
vious in Fig. [6] because of the large velocity range (36 km s~^) 
summed in this plot. Considering the overall Galactic rotation with 
roughly circular motions around the Galactic Centre, gas along the 
line of sight to the Galactic Centre will be moving transverse to 
the line of sight, with zero radial velocity. However, the motions 
within the Galaxy are not exactly circular, and the absorption fea- 
tures at diff erent velocities are associa t ed with intervening G alac- 
tic features dCreaves & Winiamslll994IWirstr6m et al.ll2010h . The 
lines are in absorption, as this intervening low density gas has lower 
excitation temperature than the background brightness temperature, 
with excitation te mperature T^x of a few K fo und for CS in sight- 
lines to Sgr B2 bv lGreaves & Winiamsl ( fl994) . 

We show in Fig.[8]the image of the peak velocity from N2H"'" 
emission data cube. This line generally shows a well defined peak 
(see also Fig. I12t and so provides a representation of the princi- 
pal kinematic velocity at each position in the CMZ. The bulk of 
the emission peaks at Vlsr ~ +50 km s^^, however there is a 
clear gradient in Galactic longitude of ~ 200 km s^^ evident, ex- 
tending from hlOO km s~^ at I = +1.5° to ~ -100 km s"^ 

at I = —0.5°. More complex kinematics is al so evident. In par- 
ticular, the twisted cold dust ring noted by Mo linari et al. I feOllh 
is apparent in the N2H"'" peak velocity image. This elliptical dust 
ring, ~ 100 x 60 pc across, has been hypothesised to trace the 
location of the system of stable X2 orbits in the barred Galactic 
potential. It is apparent in a dust column density image (see fig. 
4 Molinari et al. 201 1), and shows a coherent velocity structure a s 
traced in the emission ranee of the CS 1-0 line jTsuboi et al.l 19991) . 
The kinematics of this ring are more clearly shown in the N2H"'" 
peak velocity image shown in Fig. [8] however At / = —0.5° 
the western end of the ring is most clearly seen, with velocities 
smoothly changing from —20 to —100 km s^^ around the ring's 
extremity. The twist in the ring occurs around I ~ +0.1° as it 
crosses over itself. The eastern side, around I ~ 0.5 — 0.8°, is less 
readily identified in the velocity images, but the emission veloc- 
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ity nevertheless varies smoothly across it, from 20 — 50 km s~^ 
across the centre of the ring, to ~ +100 km s~^ at its eastern end. 
The ring twists about the Gala ctic plane, rising abov e and below 
it twice, for one complete orbit. iMolinari et al] ^201 it) propose that 
this ring is rotating about the centre of the Galaxy with an orbital 
velocity of 80 km s~^, tracing the stable, non-intersecting X2 orbits 
whose major axis is perpendicular to that of the Galactic bar, and 
processes as the the bar does. 



3.4 Integrated line emission 

For quantitative analysis of the line data cubes, it is useful to in- 
tegrate the emission over the velocity range. As the velocity range 
is large (Table [2] and subsection 13.31 above), the integrated emis- 
sion images are sensitive to noise and problems with the spectral 
baselines, as noted in subsection l3.1l This is a well known problem 
in the analysis of spe ctral-line data , as discus sed for example by 
Ivan Gorkom & EkersI ( IT989,) and .Rupenl jl999l) . 

The major problem of simply integrating the emission is spuri- 
ous stripes, in the longitude and latitude scanning directions, caused 
by spectral baseline ripples, dominated by the small fraction of data 
taken during the worst weather. A further problem is that extra 
noise is added from including parts of the spectra where there is 
no signal. 

A simple method that we adopt here is to integrate the emis- 
sion above a fixed cutoff level at which the emision is significant. 
We use a 3-a cutoff, the same used in the peak emission plots (Figs. 
[2]to[5} to define significant emission. This is quite effective in re- 
moving the spurious features in the integrated emission, but not 
totally effective as the non-gaussian statistics of the baseline errors 
means there is still a small tail to the distribution of spurious fea- 
tures above the 3 a level. Some of these artifacts can be seen in 
Figs.[2]to|5] 

However, a major problem with this simple method of inte- 
ration above a cutoff leve l is a bias in the integrated emission 
van Gorkom & Ekerdl 19891 : iRuperJI 19991) . with real emission be- 
low the cutoff being missed. Because of this, we have investigated 
two other methods of processing to obtain integrated emission. 

An alternative to the fixed cutoff level, is to integrate the emis- 
sion within a mask defined by the data smoothed spatia lly and/or 
over velocity jvan Gorkom & Ekers|[T989l : lRupenlll999l) . We have 
used as the mask a wavelet reconstruction of the cu be, implemente d 
as part of the Duchamjjf] source finding package ('Whiting| l2068h . 
This gives a model spectrum at each pixel which contains a range 
of velocity scales, but not the smallest scale, so has reduced noise, 
and so can be extended to lower brightness levels than the original 
data. This gave quite similar results in integrated emission images 
to the 3-a cutoff version, with slighly more integrated emission, 
but also slightly more artifacts. It was decided that any improve- 
ment was not worth the extra complexity of the processing. 

We have also tested different methods of processing the data 
cubes, with a spectral baseline correction, to give a filtered version, 
to help the problematic low level ripples. The baseline was deter- 
mined by an iterative process that clips the data above the 3 a level, 
replaces these clipped pixels with zero, and then smooths the data 
with a boxcar filter of width 19 pixels (around 34 km s~^). The 
baselines in parts of the data cube with significant emission (above 
the 3 a level) are unaffected. For parts of the data cube without 
significant emission, the filtering removes low level offsets to make 



the mean zero. This is similar to a high order polynomial fit to the 
spectral baselines (after the exclusion of significant line emission). 
The baseline corrections are a few times 10 mK, or around half 
the noise level in the original cubes (Table |2]l. In practice, the in- 
tegrated emission determined from the filtered cubes (after a 3-a 
cutoff, with the reduced level of a) was less than that from the un- 
filtered cubes with the original 3-a cutoff. That is, the effect of inte- 
grating to a lower level is overwhelmed by the effect of the filtering 
biasing the spectra down. 

In practice, then, we use the simple integration of emission 
above a fixed 3-a cutoff level, with the caveat that this will miss 
some flux from low-level emission. This is a tradeoff to remove the 
obvious artifiacts. However, since the baseline stability in the data 
is a limitation, such very low-level broad emission is not measured 
accurately anyway, so this a limit to the original data, rather than 
just a problem introduced by the processing. 

The integrated emission images show similar features to the 
the peak emission images Figs.|2]to|5] and look to the eye much like 
the peak emission images, but with less contrast due to overlapping 
structures, so are not plotted here. 



3.5 Principal Component Analysis 

The different spectral lines in the CMZ show similarities as well as 
differences (Figs.|2]to[5]l. It would be useful to identify and quantify 
these features of the data, in a simple objective manner. One useful 
technique to do so is Prin cipal Component Analysis (PC A), see e.g 
lHever&Schloerbl l ll997h . 

This describes the multi-dimensional data set by linear combi- 
nations of the data that describe the largest variance (the most sig- 
nificant common feature) and successively smaller variances (the 
next most significant features). 

In the context here of integrated emission images, we can 
use PCA to describe the images, with a smaller set of images 
wh ich contain the most sig nificant features. This has been used 
by lUngerechts et al] ( Il997 ) for the OMC-1 ridge, and mor e re- 
cently by us for the G333 molecular cloud jLo et al.ll2009l) and 
the Sgr B2 area with 3-mm and 7-mm molecular lines Ijones et al] 
2008c, 20 ID . 

We have implemented the PCA processing in a PYTHON 
script, with the PCA modulf^l and pyFIT^"] to read and write 
the FITS images. 

Since the PCA finds features in normalised versions of the 
different input data sets, it does not work well with low signal to 
noise data, as the noisy data is scaled up, and the PCA will be 
dominated by spurious features. So, we restrict our PCA analysis 
to the strongest eight fines (HCN, HCO+, HNC, HNCO, N2H+, 
SiO, CH3CN and HC3N) out of the total twenty lines, and an area 
around Sgr B2 and Sgr A with the strongest emission. This area was 
defined by a mask where the integrated emission of the N2H^ line 
was stronger than 10 K km s~^ . The normalisation of the integrated 
images makes versions that are mean zero and variance unity. 

The images of the four most significant components of the 
integrated emission are shown in Fig. |9] These four components 
describe respectively 82.6, 10.6, 2.8 and 1.7 percent of the variance 
in the data. These components are statistical descriptions of the in- 
tegrated line images, not necessarily physical components of the 
CMZ. However, they do highlight the physical features in a useful 



http://www.atnf.csiro.au/people/Matthew.Whiting/Duchamp/ 



http://folk.uio.no/henninri/pca_module/ 
^ ^ http ://w w w. stsci . edu/resources/soft ware. hardware/pyfits/ 
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Figure 6. The velocity structure in the CMZ as shown by the HCN line channel images, in 12 panels over 2 pages, averaged over 36 km s ^ ranges, with the 
central velocity noted in the top left of each panel. The greyscale range is to 2 K in T^. 
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Figure 7. Position-velocity diagram of the CMZ, from HCN data, sliown as peak brightness temperature, within the three latitude ranges -0.29 to -0.12 deg. 
(bottom), -0. 12 to 0.04 deg. (middle) and 0.04 to 0.21 deg. (top). Note the velocity gradient, and the absorption features at around —52, —28 and —3 km s~ ^ . 
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Figure 8. The peak velocity image, from the N2H"'" line, with the scale in km s^^, showing the overall velocity gradient in longitude, but also more complex 
kinematics. The 0.25 K contour from N2H+ peak brightness is plotted as a guide to the features. 
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way. The normalised projection of tlie molecules onto the principal 
components are shown in Fig. [TO] with successive pairs of compo- 
nents. 

The first principal component (Fig. |9]l shows emission com- 
mon to all eight strong lines, or a weighted average integrated emis- 
sion. It describes a large fraction (82.6 percent) of the variance, 
which means that all eight lines are quite similar in overall mor- 
phology. The majo r features are the ridge of emission in Sgr B2 
jjones et al. 2008a c) and the dense cores in Sgr A. 

The second principal component (Figs.l9landll0t shows the 
major difference (10.6 percent of the variance) among the eight 
lines, which is differences in the ridge of emission in Sgr B2 and 
the dense cores in Sgr A, relative to that in the first component. The 
lines of HCN, HCO+, HNC and (to a lesser extent) N2H+ show 
a relative lack of emission, which we attribute (see below, subsec- 
tion l3.6b to the brightest areas being optically thick. The other lines 
(SiO, HNCO, HC3N and CH3CN) are, relative to the first compo- 
nent, enhanced at these peaks, which is partly because they are less 
likely to be optically thick, but it is also likely that these molecules 
are preferentially found in dense cores. 

The third principal component (Figs. |9l and llOb shows small 
differences (2.8 percent of the variance) in the Sgr B2 ridge, be- 
tween most significantly SiO and HNCO compared to CH3CN. 
The fourth principal component shows small differences (1.7 per- 
cent of the variance), most significantly between HNCO and SiO. 
Although the fourth component describes only a small amount of 
the variance, it is striking that it still appears to have physical 
significance, as t he positive feature corresponds well to the north 
cloud of Sgr B2 jjones et alJboOS where HNCO is enhanced 
jMinh etal.lll998h . 

As Sgr B2 is one of the strongest features in the CMZ, and also 
a chemically complex region, some of the most significant structure 
in the PCA components here at is in the Sgr B2 region. For PCA 
analysis of t he Sgr B2 a rea with more line s at 3-mm and 7-mm see 
Ijones et ^ j2008cl) and lJones et al. U201lh respectively. 

Fig. 1 101 shows the values of the PCA component vectors, for 
the first four components, plotted as pairs (for convenience of plot- 
ting multi-dimensional vectors in two dimensions). These vectors 
describe how the normalised integrated images of different lines 
are composed of the sum of different scalar amounts (the compo- 
nents PC 1, PC 2 etc) of the principal component images (Fig.[9l(. 
All eight lines have distributions which are made up of positively 
scaled amounts (PC 1, ~ 0.3) of the first principal component, and 
further made up of positive or negative amounts (PC 2, etc) of the 
further principal components. 

We also note that the results from PCA analysis (Figs. |9] and 
llOb are robust to the potential problems of integrating the emission 
discussed above in subsection |3.4| The PCA components found are 
almost identical using integrated emission with the 3-cr cutoff here, 
to that from the integrated emission of the baselevel-filtered cubes. 

3.6 Line Intensities and Line Ratios 

The previous subsection l |3.5t on PCA shows there are significant 
differences between the morphology shown by the integrated emis- 
sion of different lines (at least for the eight strongest lines used for 
the PCA). One of the most commonly used analysis tools for spec- 
tral lines is to consider their line ratios, which we do in this sub- 
section in order to quantify the behaviour of the principal emitting 
species further across the CMZ. 

The ratio of the integrated emission of the ^^C to ^'^C isotopo- 
logues of HCN, HCO"*" and HNC lines confirms that the emission 



of the C isotopologue is generally optically thick. It is usually 
much l ess than 24, the isotopic ratio for [^^C/^'^C] at the Galactic 
Centre jLanger & Penziasll990l) . and the value that would be found 
for optically thin emission. This suggests that the differences in in- 
tegrated line distribution of HCN, HCO^ and HNC, as shown by 
the second principal component (subsection |3.5l Fig.|9l(, are related 
to optical depth effects. 

In order to present an analysis for a manageable part of this 
extensive data set we have obtained on the CMZ we have selected 
representative regions within the CMZ and examined the line emis- 
sion through averaged profiles for these apertures. The apertures 
chosen are listed in Table [3] and plotted in Fig. [TT] They include 
the entire region of the CMZ we have mapped, and the four main 
emission 'cores'; Sgr A, Sgr B2, Sgr C and G1.3. We also include 
a smaller region around the core of Sgr B2. Table[3]lists the centres 
and sizes of these apertures as well as the velocity ranges chosen 
for determining the integrated line emission and the line ratios from 
each aperture. These apertures also provide a basis for comparing 
the molecular emission from the centre of the Galaxy to that from 
the central regions of other galaxies, where generally the emission 
is spatially unresolved and often only the very brightest molecu- 
lar lines (i.e. from HCN, HNC and HCO^) can be seen. We can 
compare the equivalent parameters found for the entire CMZ, and 
with the distinct emission regions within it, to examine how they 
might vary across the central regions of the Galaxy. We can also 
determine the same parameters in the optically thin isotopologue 
lines in the CMZ (i.e. H^^CN, HN"C and H^^CO+) to examine 
whether the optical depth of particular lines may be distorting any 
interpretation being placed upon a given line ratio. 

In addition to the HCN, HNC and HCO+ lines, as listed in 
Table[T] there are other bright lines whose distribution is extended 
across the CMZ. We show the profiles for the CH3CN, HNCO, 
HC3N, N2H+ and SiO lines in Fig. [121 together with those for 
HCN, HNC and HCO+ (and their isotopologues) in Fig. [13] For 
each aperture the profiles of all lines are rather similar and are very 
broad (over 100 km s^^ wide). Some profiles show several absorp- 
tion features, for instance at —52, —28 and —3 km s~^ attributable 
to absorption by cold, foreground gas, likely in spiral arms along 
the sight l ine (e.g. as al s o seen in the AST/RO CO 4-3 datacube of 
the CMZ: ' Martin et al.l ( l2004l) ). The brightest of the lines is HCN, 
but in some instances other lines are comparably bright (e.g. HNCO 
in the Sgr B2 Core). The brightness of HCN is generally twice 
as strong as HCO"'', the second brightest of these lines, and three 
times stronger than HNC. Their isotopologue lines are significantly 
weaker, and we use the ratios with the main lines to derive optical 
depths in subsection l3.7l 

Table [4] presents the integrated line fluxes for these apertures, 
in K km s~^ and corre cted for the extended aperture efficiency of 
0.65 ( Ladd et aLll2005h . For comparison with extragalactic line lu- 
minosities L' , we integrated the brightness temperatures spatially 
and spectrally to obtain values in units of Kkm s^^ pc^ in the dif- 
ferent apertures listed in Table [5] As a reference, we also derived 
L'cn using the i ntegrated intensity map of the Columbia survey 
jPame. Hartman n. & Thaddeus Dame et al.). Over the entire CMZ, 
L'qq amounts to 2.2 x 10*^ Kkm s~^ pc^. Relative to CO, the dense 
gas tracers exhibit Ldonsc/^co luminosity ratios of 0.095 (HCN), 
0.010 (H^^CN), 0.029 (HNC), 0.002 (HN^C), 0.057 (HCO+), 
0.002 (H^^CO+), 0.005 (CH3CN), 0.014 (HNCO), 0.007 (HC3N), 
0.014 (N2H+), and 0.008 (SiO). 

Table[6]then presents the relative line ratios with respect to the 
(brightest) HCN line. The HCN, HNC and HCO+ lines may be op- 
tically thick, and so to determine whether this influences the derived 
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Figure 9. The first four principal component images derived from the eight strongest lines (HCN, HCO+, HNC, HNCO, N2H+, SiO, CH3CN and HC3N), 
which describe 82.6, 10.6, 2.8 and 1.7 percent of the variance respectively. The first principal component describes the common features of the eight lines, and 
the second principal component describes the most significant differences between the lines. The further principal components describe successively smaller 
differences. As the PCA requires good signal-to-noise data, we have restricted the lines used to the brightest and the area to that defined by the N2H+ emission 
(see text), as shown by the contour 
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Figure 10. The component vectors of the integrated Hne images in the decomposition of the data into the PCA images, shown with successive pairs of PCA 
images. HCN, HCO+ and HNC have similar vectors so the labels overlap. These vectors describe how the normalised integrated images of different lines are 
composed of the sum of different scalar amounts of the principal component images (Fig.|9). For example, all eight lines have distributions which are made 
up of positively scaled amounts (PC 1, ^ 0.3) of the first principal component, but are further made up of positive or negative amounts (PC 2) of the second 
principal component etc. 



Table 3. Apertures selected for analysis. The areas are rectangles in Galactic coordinates, centred on I and b, and using the velocity ranges indicated to include 
significant line emission. The projected area on the sky assumes a distance of 8.0 kpc. 



Source 


I 


b 


Width 


Height 


Area 


Area 


Min. vel. 


Max. vel. 




deg 


deg 


arcmin 


arcmin 


arcmin^ 


pc2 


km s~^ 


km s~^ 


CMZ 


0.545 


-0.035 


151.0 


29.8 


4 500 


24400 


-220 


220 


SgrA 


0.040 


-0.055 


30.8 


13.4 


413 


2 240 


-170 


220 


SgrB2 


0.735 


-0.045 


29.0 


18.2 


528 


2 860 


-130 


210 


SgrC 


-0.460 


-0.145 


20.0 


17.0 


340 


1 840 


-200 


170 


G1.3 


1.325 


0.015 


24.2 


20.6 


499 


2 700 


-90 


220 


Sgr B2 Core 


0.675 


-0.025 


6.2 


6.2 


38 


210 


-20 


200 



line ratios we also calculate tlie isotopologue ratios with H CN 
in Table [7] alongside the main line ratios with HCN. While the 
HNC/HCN ratio is found to be roughly the same as HN^C/H^^CN 
(~0.3) we find that the HCO+/HCN ratio (~0.6) is typically 2-3 
times as high as the H^^C0+/H"CN ratio (~ 0.2-0.3). As we will 
see below in the analysis of the optical depths, this arises because 
the HCN line is moderately optically thick (with a similar value to 
HNC) wliile the HCO"'" line has optical depth around unity. 

Molecules like CH3CN, HNCO and HC3N and N2H+ are 
generally found in dense molecular cores or cold dust cores where 
massive star formation has recently been initiated. The ratios of the 
lines from these molecules to HCN (Table |6} is seen to be more 
variable between the apertures than they are for HNC and HCO^ , 
with, in particular, Sgr B2 showing significantly larger values of 
the ratios than the other locations (typically 3-5 times larger). SiO, 
which is more sensitive to the presence of shocks, shows less vari- 
ation between the apertures, but still has a ratio twice as large for 
Sgr B2 as in Sgr A and Sgr C. However G1.3 has a similar ratio to 
Sgr B2, suggesting the influence of shocks in this region too. 

Table[8]lists the maximum and minimum line ratios found for 
HCN, HNC and HCO+, as well as for their ^^C isotopologues, 
as a function of velocity between the apertures, as opposed to the 
integrated line fluxes. The variations between these extrema can 
be seen in Fig. [14] (note that line ratios are only shown when 
the individual lines have S/N > 5 for a given velocity channel). 
Within the core of each line the ratios are generally constant, aside 



from the velocity ranges where foreground absorption is evident 
(see above). The isotopologue line ratios are generally similar for 
H"CN / HN"C compared to HCN/HNC. However for H"CN / 
H^"^CO+ the isotopologue line ratio is typically twice the value of 
HCN/HCO^, as the velocity varies. This result, for the variation 
in ratio as a function of velocity, is also consistent with behaviour 
found for the ratio of the integrated fluxes of these lines. 

However, at the highest red-shifted velocities for emission 
(> 100 km s^^), the behaviour changes, with the ratio HCN / HNC 
increasing significantly, from a value of ~ 2 over most of the pro- 
file, to ~ 10 at its high-velocity limit. This trend is also apparent 
in the "C forms, H"CN/HN^^C, though lower S/N limits the ve- 
locity range over which the comparison may be made. We note, 
however, that the HCN / HCO"'" ratio does not show any evidence 
for increasing at these high red-siiifted velocities, suggesting that 
this variation is a result of the relative strength of the HNC line de- 
creasing with increasing velocity, while the relative intensity of the 
HCN and HCO+ lines remains constant. 



3.7 Optical Depth Analysis 

Where the H^CN, HN^^C and H^^CO+ lines are well detected 
(S/N > 5) we may use the ratio with the main isotopologue to 
determine the optical depth of the emission, as a function of ve- 
locity across the line ratio. We list in Table |9] the maximum and 
minimum values of the ratios found for the isotopologue pairs (i.e. 
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Table 4. Mean integrated line fluxes (K km s^^), as jTMsdV (i.e. corrected for aperture efficiency) averaged over the apertures and integrated over the 
velocity ranges given in Table[3] For each molecule the first line gives the integrated flux and the second line its statistical 1 cr error Note that for the weaker 
lines the uncertainty in where to place the baseline across the wide line profi l e means tha t the formal S/N for the integrated flux is often low, even though the 
lines themselves are clearly detected. The CO fluxes are obtained from i200ll jPame et al.h . 



Source 


CO 


HCN 


H^^CN 


HNC 


HN^^C 


HCO+ 


h1^CO+ 


CHgCN 


HNCO 


HC3N 


N2H+ 


SiO 


CMZ 


914 


87 


9 


26 


2 


52 


2 


4 


13 


6 


13 


7 






1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


SgrA 


1200 


140 


18 


48 


4 


85 


4 


8 


10 


10 


25 


9 






1 


2 


2 


4 


2 


3 


2 


4 


2 


2 


4 


SgrB2 


1 170 


142 


19 


45 


5 


86 


5 


11 


40 


18 


33 


17 






1 


1 


2 


2 


1 


1 


1 


2 


1 


1 


2 


SgrC 


1060 


100 


9 


35 


3 


61 


1 


5 


5 


6 


13 


3 






2 


1 


2 


3 


2 


3 


1 


3 


2 


1 


4 


G1.3 


1270 


132 


12 


30 


3 


73 


2 


4 


20 


8 


17 


17 






1 


1 


1 


1 


1 


2 


1 


2 


1 


1 


2 


Sgr B2 Core 


1610 


199 


31 


76 


11 


123 


10 


33 


101 


48 


68 


34 






3 


3 


5 


3 


2 


4 


3 


3 


3 


2 


2 



Table 5. Line luminosities (x 10^ K km s^^ pc^), / TMBdVdA integrated over each aperture and the velocity ranges given in Table[3] For each molecule 
the first fine gives the line luminosity and the second line its statistical 1 a error. Note that for the some lines the uncertainty in where to place the baseline 
across the wide li ne profile means th at the formal S/N for the luminosity is poor, even though the fines themselves are clearly detected. The CO luminosities 
are obtained from I2OOII jPame et alj) . 



Source 


CO 


HCN 


H^^CN 


HNC 


HN^-^C 


HCO+ 


Hi3CO+ 


CH3CN 


HNCO 


HC3N 


N2H+ 


SiO 


CMZ 


2 230 


212 


22 


64 


5 


128 


5 


11 


32 


16 


32 


17 






1 


1 


2 


2 


1 


2 


2 


2 


1 


2 


3 


SgrA 


268 


31.4 


4.1 


10.8 


0.8 


19.1 


0.9 


1.7 


2.3 


2.1 


5.7 


2.1 






0.2 


0.4 


0.5 


0.9 


0.5 


0.8 


0.4 


0.8 


0.3 


0.5 


0.8 


SgrB2 


335 


40.6 


5.4 


12.9 


1.5 


24.6 


1.5 


3.1 


11.4 


5.2 


9.3 


4.8 






0.3 


0.3 


0.5 


0.5 


0.3 


0.3 


0.3 


0.5 


0.4 


0.2 


0.4 


SgrC 


195 


18.4 


1.7 


6.4 


0.6 


11.3 


0.2 


0.9 


1.0 


1.1 


2.4 


0.6 






0.3 


0.2 


0.3 


0.5 


0.4 


0.5 


0.3 


0.6 


0.3 


0.2 


0.8 


G1.3 


343 


35.6 


3.4 


8.0 


0.7 


19.6 


0.5 


1.2 


5.3 


2.2 


4.5 


4.7 






0.4 


0.2 


0.3 


0.4 


0.2 


0.5 


0.2 


0.5 


0.3 


0.3 


0.5 


Sgr B2 Core 


33.4 


4.14 


0.64 


1.59 


0.24 


2.56 


0.20 


0.68 


2.09 


1.00 


1.42 


0.71 






0.06 


0.06 


0.11 


0.06 


0.05 


0.08 


0.06 


0.06 


0.06 


0.05 


0.05 



Table 6. Ratios of integrated fines fluxes with HCN. For each molecule the first line gives the ratio of that line and the HCN 1-0 flux and the second line the 
1 a en'or. For HCN itself the errors indicate the flux eiTor range for that molecule. 



Source HCN H^'^CN HNC HN^'^C HCO+ H1^C0+ CH3CN HNCO HC3N N2H+ SiO 



CMZ 


1.00 


0.10 


0.30 


0.03 


0.60 


0.03 


0.05 


0.15 


0.07 


0. 


15 


0.08 




0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0. 


01 


0.02 


SgrA 


1.00 


0.13 


0.34 


0.03 


0.61 


0.03 


0.05 


0.08 


0.07 


0. 


18 


0.07 




0.01 


0.01 


0.02 


0.03 


0.02 


0.02 


0.01 


0.03 


0.01 


0. 


02 


0.03 


SgrB2 


1.00 


0.13 


0.32 


0.04 


0.61 


0.04 


0.08 


0.28 


0.13 


0. 


23 


0.12 




0.02 


0.01 


0.02 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0. 


01 


0.01 


SgrC 


1.00 


0.09 


0.35 


0.03 


0.61 


0.01 


0.05 


0.05 


0.06 


0. 


13 


0.03 




0.04 


0.01 


0.02 


0.03 


0.03 


0.03 


0.01 


0.03 


0.02 


0. 


02 


0.04 


G1.3 


1.00 


0.09 


0.22 


0.02 


0.55 


0.01 


0.03 


0.15 


0.06 


0. 


13 


0.13 




0.02 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.02 


0.01 


0. 


01 


0.02 


Sgr B2 Core 


1.00 


0.16 


0.38 


0.06 


0.62 


0.05 


0.16 


0.51 


0.24 


0. 


34 


0.17 




0.03 


0.02 


0.03 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0. 


02 


0.01 
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1.5 



0.5 

GLON (degrees) 



359.5 



Figure 11. The six areas selected for analysis: ai'ound Sgr A, Sgr B2, Sgr C, G1.3, a smaller area at the Sgr B2 core, and the larger CMZ ai'ea, plotted on the 
peak HCN emission. 



Table 7. Ratios of integrated line fluxes for HNC and HCO+ with HCN, together with the same ratios measured for their isotopologues. 



Source 


HNC / HCN 


HNi3C/H"CN 


HCO+ / HCN 


Hi3C0+ /Hi^CN 


CMZ 


0.30 ± 0.01 


0.25 ±0.10 


0.60 ± 0.01 


0.25 ± 0.09 


Sgr A 


0.34 ± 0.02 


0.20 ± 0.24 


0.61 ± 0.02 


0.23 ± 0.22 


SgrB2 


0.32 ± 0.02 


0.28 ±0.11 


0.61 ±0.01 


0.27 ± 0.08 


SgrC 


0.35 ± 0.02 


0.35 ± 0.34 


0.61 ± 0.03 


0.11 ±0.29 


G1.3 


0.22 ± 0.01 


0.22 ±0.12 


0.55 ± 0.01 


0.15 ±0.17 


Sgr B2 Core 


0.38 ± 0.03 


0.37 ±0.13 


0.62 ± 0.02 


0.31 ±0.15 



Table 8. Maximum and minimum values of selected line ratios with the velocities (V^sr) which these extrema occur. Note that line intensities need to be 
> 5(7 per channel to be used for this analysis. No data meets this criteria for HN^^C in Sgr A. 



Source 




HCN / HNC 






h"cn/hn '^'^c 






HCN / HCO+ 






H^^CN/H 


13CO+ 






Max. 


Vel. 


Min. 


Vel. 


Max. 


Vel. 


Min. 


Vel. 


Max. 


Vel. 


Min. 


Vel. 


Max. 


Vel. 


Min. 


Vel. 




ratio 


kms~^ 


ratio 


kms~^ 


ratio 


kms~^ 


ratio 


kms~^ 


ratio 


kms^^ 


ratio 


kms^^ 


ratio 


kms^^ 


ratio 


kms^^ 


CMZ 


7.5 


183 


1.5 


-49 


8.2 


77 


0.7 


-116 


2.5 


-192 


0.9 


223 


6.6 


101 


2.6 


52 


Sgr A 


7.8 


97 


0.7 


-1 










5.1 


-7 


0.9 


-47 


3.8 


68 


2.6 


45 


SgrB2 


9.9 


127 


1.2 


-41 


5.3 


94 


2.2 


54 


2.6 


-94 


0. 6 


-36 


6.3 


4 


2.0 


54 


SgrC 


6.8 


81 


1.1 


-50 


2.5 


-39 


0.4 


-92 


3.0 


-54 


0.7 


-49 


2.5 


-56 


2.3 


-59 


G1.3 


15.9 


139 


1.3 


-34 


6.2 


101 


1.4 


34 


2.1 


75 


1.0 


212 


5.4 


94 


3.7 


70 


Sgr B2 Core 


5.6 


117 


1.4 


50 


5.5 


86 


1.3 


54 


3.0 


-1 


1. 3 


143 


4.7 


86 


1.4 


52 



Table 9. Maximum and minimum values of the ratios of the line intensities, per velocity channel, for HCN, HNC and HCO+ with their '^C isotopologues, 
together with the velocities (Vlsr) ^t which these extrema occur. Note that all line intensities need to be > 5cr per channel to be used for this analysis. No 
data meets this criteria for HN^'^C in Sgr A. 



Source 




HCN/H 


i^CN 






HNC/HN13(;. 






HCO- 








Max. 


Vel. 


Min. 


Vel. 


Max. 


Vel. 


Min. 


Vel. 


Max. 


Vel. 


Min. 


Vel. 




ratio 


km s~^ 


ratio 


km s^^ 


ratio 


km s^^ 


ratio 


km s"'^ 


ratio 


km s~ 


^ ratio 


km s~^ 


CMZ 


22.3 


139 


4.6 


-1 


23.3 


77 


2.1 


-114 


44.1 


101 


1 0.1 


-28 


Sgr A 


11.1 


-37 


1.5 


-1 










20.3 


68 


1 2.4 


46 


SgrB2 


15.8 


127 


2.6 


-47 


14.1 


94 


6.2 


52 


38.6 


103 


8.4 


52 


SgrC 


29.9 


-78 


2.9 


-54 


8.9 


-39 


2.2 


-83 


14.2 


-59 


7.0 


-56 


G1.3 


18.9 


143 


3.3 


-36 


16.4 


81 


3.1 


26 


27.4 


94 


1 6.6 


83 


Sgr B2 Core 


10.4 


110 


4.1 


46 


13.9 


83 


4.4 


54 


22.3 


86 


4.0 


52 
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Figure 12. Profiles of five molecular lines in six selected apertures in the CMZ, as defined in Table|3]and Fig. 1111 From left to right, the lines are (see Table[T) 
from the CH3CN, HNCO, HC3N, N2H+ and SiO molecules. The apertures run from top to bottom and are: the entire CMZ, Sgr A, Sgr B2, Sgr C, G1.3 and 
the Sgr B2 Core. Intensities are the averaged pixel values in each aperture, in Kelvin, and are (i.e. without beam efficiency correction). The x-axis is the 
velocity, Vj^sR^ running from -200 to +200 km s~^. 
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Figure 13. Profiles of the three bright molecular transitions, the HCN, HNC and HCO+ J=l-0 lines in the Central Molecular Zone (left, middle and right 
columns, respectively). Overplotted as dashed lines are the corresponding isotopologues; H'^^CN, HN^^C and H^''CO+ J=l-0, multiplied by a factor of 3 for 
clarity. From top to bottom are the six selected apertures (see Table[3]and Fig. II It : the entire CMZ, Sgr A, Sgr B2, Sgr C, G1.3 and Sgr B2 Core. Intensities 
(TJ) are the averaged values over each aperture and velocities are Vi^sB. and run from -250 to +250 km s~^. 
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Figure 14. Line ratios, as a function of velocity, for HCN/HNC (left) and HCN/HCO+ (right), for the six selected apertures in Table[3](the entire CMZ, Sgr A, 
Sgr B2, Sgr C, G1.3 and the Sgr B2 Core). The main line is shown in soHd and the l^C-isotopologue ratios (Hi^CN/HN^C and H13cN/H"C0+) are 
dashed. Line ratios are only shown when the intensities of both lines are > 5 it in the relevant velocity channel. Note that for HN^^C in Sgr A this results in 
no useable isotopologue data (and just 2 data points for H^^CO+ in Sgr C). The x-axis is the V^sR velocity, running from -150 to +150 km s~^. Note that 
the y-axis ranges are the same in all plots except for Sgr B2 and G1.3 (left column). 
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HCN/H^^'CN, HNC/HN"C and nCO+IW'^CO+) across their line 
profiles, together with the velocities at which these extrema occur. 
Again, these numbers are filtered to require the S/N in each velocity 
channel to be > 5cr. The variation of each line ratio with velocity 
can be followed in Fig. [15] The ratios are generally well less than 
24, the abun dance ratio for ^^C to ^^C determined for the centre of 
the Galaxy jLanger & Penziajl 19901) , indicating that the ^^C iso- 
topologues of these species must be optically thick. However, the 
isotopologue ratios are also significantly greater than unity, which 
also implies that the ^''C lines cannot be strongly optically thick. 
We make use of this behaviour to provide an approximate solution 
of the radiation transfer equation below. We apply a standard anal- 
ysis technique, which we summarise here. 

The radiative transfer equation has the general solution for the 
intensity of a velocit y channel in a spectra l line known as the de- 
tection equation (e.g. istahler & Pallal20o3) : 

Ta = T\lri = j{J^(T^^) - .h{TBG)\(\ - exp(-r,.)) 

where is th e measured inten sity, 77 the extended beam efficiency 
(taken as 0.65, iLadd et al.l20o3) . / is the beam filling factor for the 
emission, Ji/(T) = [/i://fc]/[exp(/i://fcT) — 1], with Tea, being 
the excitation temperature and Tbg being the temperature of the 
CMBR (i.e. 2.726 K). r„ is the optical depth of the emission. 

For a line isotopologue pair, say HCN and H^"^CN, we as- 
sume the former is optically thick, but the latter is optically thin, 
as inferred above. We also assume they have the same excitation 
temperature and beam filling factor. Furthermore, 

the abundance ratio of the ^^C isotope to ^''C, which we have taken 
to be 24 (Langer & Penzias 1990), assuming no isotope fractiona- 
tion in the molecules. 

Then we see, from application of the detection equation to 
HCN and H"CN, that rj^iacjv = TX{H^^CN)/TX{HCN), 
with thcn then determined by the preceding formula. 

We hence apply this analysis to determine the optical depth 
of the HCN, HNC and HCO"'" lines, as a function of velocity. Fur- 
thermore, from stan dard molecular radiative transfer theory (e.g. 
lOoldsmith & Langer 1999) we may show that the column density 
of the upper level of each transition is given by 

Nu = [^-Kk/AhWu^ /c^]Ta[t/1 - exp(-r)](5V 

where k, h and c are the well-known physical constants, v is the 
frequency of the transition, A its radiative decay rate (Einstein 
coefficient) and 5V the channel velocity spacing. The optically 
thin case (r <C 1) simply has the optical depth correction factor 
r/(l — exp(— r)) set to unity. 

Once A'^^u has been determined the total molec- 
ular column density A'^ can be found by applying 
= {NuQ{T)/gu) exp{E/kn^) with ff„ = 2J 1 = 3, the 
partition function Q{T) = 2kT/hv, E the energy of the upper 
level and an assumed excitation temperature, Tex- We take the 
later as 24 K, base d on the dust temperature of the Sgr B2 envelope 
jjones et al]|201lh . Nagavama et al. (2009) determined that Tex 
has the range 20 to 35 K in the CMZ from LVG modelling of 
the CO lines. To check the sensitivity of this assumption we also 
repeated our calculations by halving and doubling the temperature. 
For the former (i.e. T = 12 K) the column densities (and masses) 
decrease by a factor ~ 0.6. For the latter (i.e. T = 48 K) they 
increase by a factor ^ 1.8. While undoubtedly variations in the 
excitation temperature across the CMZ could result in relative 
variations in these calculated quantities by factors of a few, the 



absolute determinations for the column density in each aperture 
are not likely to be in error by more than a factor 2 due to this 
assumption of a constant temperature T = 24 K 

Similarly the molecule mass M can then be determined from 
the column density TV knowing the aperture size (Table[3l(, molec- 
ular masses rrimoi, and the source distance (D = 8.0 kpc) as 
M = Nnimoi^D^. 

In TablellOlwe list the derived optical depths for each aperture 
for HCN, HNC and HCO+. We list both the maximum value of 
the optical depth found at each position (together with the relevant 
emission velocity) and the weighted mean optical depth. This is de- 
termined from the value of the correction factor r/(l — exp(— r)) 
needed to convert the integrated line flux from that for an optically 
thin line to the measured integrated line flux after applying the op- 
tical depth correction at each velocity channel. This number thus 
provides a representative value of the optical depth for each line. 
For HCN and HNC this is typically 2-3, whereas for HCO^ it is 
around unity. On the other hand, the maximum values of the opti- 
cal depth for each line are several times their mean values. These 
are illustrated graphically in Fig. [TS] This shows, for each source 
and line, both the measured and the optical-depth corrected line 
profile. Also shown, with the right-axis as the numerical value, are 
the measured value of the isotopologue line ratio (i.e. HCN/H^'^CN 
etc.) and the value determined for the optical depth at each veloc- 
ity. The optical depth is found to be fairly constant over each line 
though note that it is higher where the absorption features are seen; 
however by applying the optical depth correction these features are 
removed from the line profiles, as can be seen by their smoothly 
passing across the affected velocities in Fig. 1151 

While the isotopologue line ratios are not found to vary much 
over the cores or blue-wings of each line, this is not the case for 
the red-shifted gas, however. For several apertures (especially the 
entire CMZ, Sgr B2 and G1.3; see Fig.|15t these ratios rise steeply 
for velocities > +50 km s~^, from values of ~ 5 to reach 20^0 
for the most red-shifted emission (see also Table [9l note that this 
analysis has only been applied to data with S/N > 5 per channel 
so this is not an artifact of low S/N at the profile edges). This is 
as expected if the optical depth is falling with increasing velocity, 
with the line ratio tending towards the optically thin value given by 
the [^^C/^'^C] isotope ratio. However, we note that the ratios deter- 
mined for the entire CMZ and Sgr B2 for the HCO^ isotopologues 
at the highest velocities are > 30, greater than the assumed iso- 
tope abundance ratio of 24 for the CMZ. We discuss this further in 
subsection |4.3l 

Finally Table [TT] presents the derived column density and 
masses for each aperture and molecule. Presented here are the up- 
per level column density, without any optical depth correction, and 
the total molecule column density, correcting for optical depth and 
assuming an excitation temperature of 24 K. The molecule masses 
are then derived for each aperture from the total column density for 
that molecule. The column densities for each molecule vary little 
with aperture, though of course the masses do as they are summed 
over the aperture. Around 50 M© of [HCN-1-HNC-1-HCO+] exists 
in the CMZ, 60 percent of this being HCN and 20 percent each 
HNC and HCO^. Furthermore, 10-15 percent of the total mass of 
the CMZ is associated with each of the regions Sgr A, Sgr B2 and 
G1.3. 

We end this section by noting that we have considered cal- 
culating the molecule abundances by making use of the sub- 
millimetre dust data to provide total column densities in each aper- 
ture, and then ratioing these values with the relevant molecular col- 
umn densities. Three such data sets have been published, at 1.1 mm 
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Figure 15. Optical depth-corrected line intensities, together with the optical depth and isotopologue ratios, as a function of velocity. From left to right are the 
HCN, HNC and HCO+ measured line protiles (soUd Une) and optically-depth corrected profile (dashed line). The intensity, T^, in Kelvin is shown on the 
left-hand axis, and the x-axis shows the velocity in Vlsr from —150 to +150 km s~^. The symbol X shows the measured isotopologue line ratio for each 
molecule (i.e. HCN/Hl^CN, HNC/HNi^C and HC0+/Hl3C0+) and the + sign shows the calculated optical depth for the HCN. HNC and HCO+ J=l-0 
lines, respectively (see text). The right-hand axis indicates the numerical value for the line ratio and optical depth. The six selected apertures (see Table|3j are, 
from top to bottom, the entire CMZ, Sgr A, Sgr B2, Sgr C, G1.3 and the Sgr B2 Core. Note that optical depths and line ratios are only calculated when the 
con'esponding lines are > 5 ct (i.e. HCN and H^'^CN > 5 o", etc.). For some lines/sources (e.g. HN^'^C in Sgr A) there are few, or no, data for which this 
condition is met. Also, for clarity, every three data points in the x-direction are binned together to provide a single point on the plots. 
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Table 10. Optical depths, for HCN, HNC, and HCO"'", derived from the isotopologue ratios assuming [^^C/^''C] = 24, given as the maximum and the weighted 
mean. The velocity listed is that at which the maximum optical depth occurs. See the text for the definition of the weighted mean value of the optical depth. 
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HCN 


HCN 


HCN 


HNC 


HNC 


HNC 


HCO+ 


HCO+ 


HCO+ 




Max. 


Vel. 


Mean 


Max. 


Vel. 


Mean 


Max. 


Vel. 


Mea n 
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km s~^ 
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T 


km s~^ 


r 


T 


km s^^ 


T 


CMZ 


5.2 


-1.1 


2.4 


11.7 


-114.1 


1.0 


2.4 


-28.4 


0.9 


SgrA 


16.5 


-1.1 


3.0 








1.9 


46.3 


1.0 


SgrB2 


9.1 


-46.6 


3.2 


3.9 


51.8 


2.5 


2.9 


51.8 


1.4 


SgrC 


8.2 
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2.2 


Il.O 


-83.1 


2.9 


3.4 


-55.8 


0.1 


G1.3 


7.4 


-35.7 


2.2 


7.6 


26.3 


1.6 


1.4 


82.8 


0.4 


Sgr B2 Core 


5.9 


46.3 


3.6 


5.5 


53.6 


2.6 


5.9 


51.8 
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(BOLOCAM. iBa Iv et al.1 IZOlOh . at 450 and SSO A^m (SCUBA 
Pierce-Price et al.ll200(]|) and at 870/im (LABOCA. ISchuUer et al] 



20091) . However close examination of these data sets reveals i 



with the determination of the baseline in an extended emission re- 
gion such as the CMZ, given the data acquisition methods applied. 
Significant levels of low-level emission are missing in the final im- 
ages, so yielding incorrect dust column determinations, even be- 
fore other uncertainties such as dust temperature and opacities are 
considered. We have therefore not attempted to determine molec- 
ular abundances from our data set, which await further work on 
the analysis of the extended dust emission across the CMZ before 
reliable values can be calculated. 



4 DISCUSSION 

4.1 CMZ Line Luminosities Compared to External Galaxies 

iGao & SolomonI ( |2004|) surveyed 52 galaxies in HCN(l-O) that 
were selected to cover a large range of FIR luminosities. Their 
sample galaxies show L'qq values in the 5 x 10® to 1.5 x 
lO^^Kkms"^ pc^ range. This is considerably more than what is 
measured in the CMZ. Even galaxies that are considered Milky 
Way equivalents, such as IC342, show over an order of mag- 
nitude larger L'qq luminosities. The reason for the difference is 
likely the edge-on orientation of the CMZ, the sample selection 
that picks galaxies with rather large FIR fluxes, as well as the fact 
that the physical size of their single dish beam projects to kpc 
size scales at the distances of the sample galaxies - much larger 
than t he ^350 pc of ou r CMZ image. The ratio of L'^q-^/L'qq 
in the lOao & SolomonI 12004) sample, however, is in the range 
0.03 - 0.18, which is very similar to what we find in the CMZ 
(0.095) . This is also in a greement with the more recent observa- 
tion of iKrips et alj ( l2008h . They observed a sample of 12 nearby, 
starburst and active galaxies and derive I/hcn/^co luminosity 
ratios of 0.1-0.5 for the bulk of their sample of 12 galaxies. The 
CMZ is at the lower end of that range and may indicate that the star 
formation activity correlates with th e availability of dense gas rel- 
ative t o CO, a trend establish ed by [Solomon, Downes & Radford 



(Il992[) and |G ao & Solomon (2004|) a nd interpret e d, e.g ., by 
iKrumholz. McKee & Tum linson (20091). iKrips et alj ( l2008h also 
measure HCO"*" and they derive 1/^00+ /-^hcn ratio in the 0.5 to 
1.6 range. Our value for the CMZ is with 0.6, at the lower end of 
that distribution. 



4.2 Interpreting Line Ratio Variations for HCN, HNC and 
HCO+ in External Galaxies 

The ratios of the fluxes of the HCO+, HCN and HNC lines are of- 
ten measured in studies of extra-galactic sources, where they may 
provide the only readily accessible probe of the dense molecu- 
lar environment ac ross the entire central regions of a Galaxy (e.g. 
iLoenen et alj2008h . In particular, the ratios of these lines have been 
used to distinguish X-ray-dominated regions (XDRs) from pho- 
todissociation regions (PDRs) (Meijerink, Spaans, & Israel 200"^, 
and applied to starburst galaxies l lBaan et alj2008l) . 

The spatially resolved ratios for these lines that we have deter- 
mined from the CMZ are therefore useful in providing a Galactic 
analogue to compare such results with. Within the CMZ it is possi- 
ble to study how the ratio may change with position, due to chang- 
ing excitation conditions, as well as to see how it appears when 
blended together in a single pixel as in an external galactic nucleus. 
We may also examine the effects of opacity on the diagnostic line 
ratios by comparing them to the ratio obtained using their optically 
thin isotopologues. The isotopologues may not be measurable in 
an external galaxy due to their weakness, yet excitation conditions 
may be also mis-diagnosed if the measured line ratios do not reflect 
the intrinsic line ratio due to one line of a pair being more affected 
by optical depth than the other. 

In such analyses of external galaxies generally the logarithm 
of the line ratios (e.g. log[I(HCO+)/I(HCN)] etc.) is presented. To 
aid this comparison we have therefore re-calculated the line ratios 
in Table|7]into their logarithm values and present these in Tablell2l 
together with the errors. For completeness we have also included 
HCO+/HNC in this table although it can be calculated from the 
other two ratios. The value for the ^''C line ratios in Table [T2l can 
be interpreted as providing the optically thin ratio for the corre- 
sponding ^^C line ratios across the CMZ. 

For log(HNC/HCN) the line ratio varies from -0.4 to -0.6 
between the apertures. The ^^C line ratios are consistent with those 
for the corresponding ^^C value, within the errors of measurement. 
However, this is not the case for log(HCO^/HCN). This varies lit- 
tle between the 6 apertures we have chosen for the ^^C line ra- 
tios, from —0.21 to —0.26. Yet for the ^'^C isotopologues the cor- 
responding value is typically a factor 2 higher (i.e. ranging from 
-0.5 to -0.6). For log(HCO+/HNC) there is a factor of 2 vari- 
ation measured between the apertures, from 0.2 to 0.4. However, 
while the measured ^^C isotopologue ratios for these species may 
appear to differ from this in Table [72l the relevant lines are too 
weak to discern any significant difference from the corresponding 
^^C ratios. To summarise, the empirical result from the CMZ data 
presented here is that the HCO^/HCN ratio needs to be determined 
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Table 11. Column density and mass for HCN, HNC and HCO^". The column densities are the mean over each aperture, and calculated as Nj—i, in the J = 1 
level without the coirection for optical depth and, as A'^ the sum over all levels, with the correction for optical depth (t/(1 — exp(— r))), applied per velocity 
channel and assuming T^^ = 24 K. The mass, M, is the sum over the aperture defined in Table[3] calculated including the correction for optical depth. Note 
if Tex is halved A'^ and M decrease by 0.6 times, while if Tex is doubled they increase by 1.8 times. 



Source 




N 


M 




N 


M 




N 


M 




HCN 


HCN 


HCN 


HNC 


HNC 


HNC 


HCO+ 


HCO+ 


HCO+ 




no corr. 


corr. 


corr. 


no corr. 


corr. 


corr. 


no corr. 


COIT. 


corr. 






1014 cm-2 


M0 


1014 cm-2 


1014 cm-2 


M© 


1014 cm-2 


1014 


M© 


CMZ 


1.5 


17 


29 


1.0 


7.0 


11 


0.74 


5.0 


8.8 


SgrA 


1.6 


23 


3.5 


0.96 


4.2 


0.64 


0.72 


5.1 


0.81 


Sgr B2 


1.7 


25 


4.8 


0.84 


9.9 


1.9 


0.74 


6.3 


1.3 


SgrC 


1.4 


15 


1.9 


0.90 


12 


1.5 


0.60 


2.8 


0.37 


G1.3 


1.4 


16 


2.9 


0.61 


5.4 


0.98 


0.65 


3.5 


0.68 


Sgr B2 Core 


1.8 


31 


0.43 


0.84 


11 


0.15 


0.76 


7.9 


0.12 



from measurements of the C isotopologues, but the HNC/HCN 
ratio can be determined using the main ^^C lines. We cannot deter- 
mine whether the I'^C isotopologues are necessary to measure the 
true HCO^/HNC ratios. However the analysis given earlier, show- 
ing that the HCO"'" optical depth is generally less than that for HNC, 
sug gests that this i s likely to be the case. 

iLoenen etal] j2008) examine the HCN, HNC and HCO+ line 
ratios in a number of IR-luminous galaxies and present plots of 
their variation for the three i2C-species ratios considered in Ta- 
blell2l They discriminate in these plots between excitation regimes 
dominated by PDRs from XDRs; HNC and HCO+ are found to be 
relatively stronger, compared to HCN, in XDRs than PDRs. They 
find that nearly all the sources they studied lie clearly in the PDR- 
regime; indeed they need to invoke an additional heating source 
("mechanical heating") even to explain the ratios under the PDR 
regime. 

Applying the results from our CMZ study actually exacerbates 
this situation. All points for the CMZ also lie clearl y in the PDR 
regime in these plots (see Fig. 1 of lLoenen et al.l20o3) , which is not 
unexpected given the W^Lq IR luminosity of the CMZ is domi- 
nated by re-processed dust emission. The CMZ data points lie in the 
same region of the phase space in these diagrams as the majority of 
the galaxies in the study. However, applying the H^^C0"''/H1^CN 
ratio as a diagnostic, instead of HCO"'"/HCN, moves the CMZ 
points even further from the dividing line between the PDR and 
XDR regimes. By analogy, this suggests that this is also likely to 
be the case for the majority of galaxies where these diagnostic ra- 
tios have been measured. 



4.3 Isotopologue Line Ratio Variations with Velocity 

As can be seen in Fig. [TS] and was discussed in subsection 13.71 the 
line ratios for the three isotopologue pairs increase significantly 
with red-shifted velocity in several locations (CMZ, Sgr B2, G1.3). 
The trend is particularly noticeable for HCO^ and HCN, and less 
so for HNC. This can be caused by decreasing optical depth of the 
12c species with increasing velocity, with the line ratio then tend- 
ing towards the [i^C/^'^C] isotope ratio (at least in the absence of 
fractionation effects). However we note that there is no trend of in- 
creasing ratio (or decreasing optical depth) seen to the blue-shifted 
side of the line profiles. 

Of the three line ratio pairs, the HC0+/H1^C0+ ratio is gen- 
erally the highest. This is also consistent with the average optical 
depth of the HCO+ line being lower than HCN or HNC (see Ta- 



blellOll. HCO"'" provides, therefore, the best probe of the optically 
thin ratio; i.e. of the [i^C/^'^C] isotope ratio. 

For the Sgr B2 and integrated CMZ apertures the 
HC0+/H1^C0+ ratio rises as high as ~ 40 at the high-velocity 
edge of the profile, considerably in excess of the [i^C/^^C] aver- 
age isotope ratio of 24 for the centre of the Galaxy. Taken by itself, 
the rising ratio with decreasing HCO^ intensity in the red-shifted 
wing could be a result of a systematic error in the baseline subtrac- 
tion of the I'^C species, making the calculated flux relatively fainter 
compared to the I'^C species than it should be. However a simi- 
lar result has also been reported by Riquel me et alj (|2pi0b) in the 
Galactic Centre from observations using the IRAM 30m telescope. 
They find similarly high HC0+/H"C0+ ratios (~ 40) at several 
selected positions in the central 5° of the G alaxy, with the ratio also 
highest for the most red-shifted gas. Riqu elme et alj (2010b) also 
measure lower values for the HCN and HNC isotopologue ratios, 
similar to our results. 

If these inferred isotope ratios of > 24 are correct then this 
red-shifted gas has undergone a different degree of nuclear pro- 
cessing that the bulk of the CMZ gas, which has a [i2c/i^C] iso- 
tope ratio of 24. i^C is primarily produced in first generation, 
massive stars on a rapid timescale. I'^C arises from CNO pro- 
cessing of i^C in lower mass stars , on a slower time-scale (e.g. 
IPrantzos. Aubert. & Audouz3 Il996h . There is also a gradient in 
[i2c/^^C] observed with distance from the centre of the Galaxy, 
rising from ^ 24 there to reach 80-90 in the solar neighbourhood 
( IWilsoiJl99Si) . lRiauelme et"^j2010bl) interpret these elevated iso- 
tope ratios resulting from infall of material to the centre of the 
Galaxy. Positions observed associated with xi orbits (which are 
parallel to the Bar) had higher inferred values of [^^C/^'^C] than 
those positions associated wi th X2 orbits (whic h are inside and or- 
thogonal to the xi orbits; see IStark et al.l jBoj) ). 

There is limited o verlap between the positions observed by 
iRiguelme et al. I bo 1 Obi) and ourselves; moreover the large aper- 
tures we have used for our analysis are in contrast to the single 
beams used by Riquelme et al. (2010b). Nevertheless, the Sgr B2 
aperture we applied includes b oth a position a t high-v elocity as- 
sociated with an xi orbit where iRiquelme et alj l l2010b l) measured 
[i2c/^^C]~ 56 and the Sgr B2 core where they measured ~ 4. 
These are similar to the ratios we measure; in the extended Sgr B2 
aperture [i^C/^^C]is ~ 40 at Vlsr = +100 km s^i, and in the 
Sgr B2 core is ~ 4 at Vlsr = +50 km s"i. The high i2c/^'^C ra- 
tio we measure at the most red-shifted velocities is thus consistent 
with the CMZ being fed with gas with [i^C/^^C] ratios typical of 
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Table 12. Logarithm of line ratios, for the integrated flux through each aperture, for HCN, HNC and HCO^, together with same ratios measured for the C 
isotopologue. Statistical errors are also listed for each ratio. 



Source 


log(HCO+/HCN) 


log(H"cO+/Hi3cN) 


log(HNC/HCN) 


log(HN"c/Hi3CN) 


log(HCO+/HNC) 


log(Hi3cO+/HNi3C) 


CMZ 


-0.22 ±0.01 


-0.59 ±0.13 


-0.52 ±0.02 


-0.61 ±0.16 


+0.30 ±0.02 


+0.02 + 0.21 


SgrA 


-0.22 ±0.02 


-0.64 ±0.29 


-0.46 ± 0.03 


-0.70 ±0.35 


+0.25 ±0.04 


+0.06 ± 0.47 


SgrB2 


-0.22 ±0.01 


-0.57 ±0.11 


-0.50 ±0.03 


-0.54 ±0.14 


+0.28 + 0.03 


-0.02 + 0.20 


SgrC 


-0.21 ±0.04 


-0.97 ±0.57 


-0.46 ± 0.04 


-0.45 ±0.30 


+0.25 + 0.05 


-0.52 ±0.65 


G1.3 


-0.26 ±0.02 


-0.83 ±0.33 


-0.65 ± 0.03 


-0.67 ±0.21 


+0.39 + 0.03 


-0.15 ±0.42 


Sgr B2 Core 


-0.21 ±0.03 


-0.50 ±0.18 


-0.42 ±0.05 


-0.43 ±0.14 


+0.21 + 0.05 


-0.07 ±0.22 



those found at larger distances from the Galactic Centre. A future 
investigation will be to examine the isotopologue line ratios at po- 
sitions across the CMZ where less nuclear-processed material may 
be falling into the Galactic Centre onto the xi and X2 orbits, so 
as to determine whether the above result is spread throughout the 
region, or if it is simply a special example. 

It would also be of interest to examine the CO isotopologue 
(i.e. ^^C^O, "C^'^O and ^^C^^^O) ratios as a function of veloc- 
ity to determine whether this behaviour is seen for these lines in 
[^^C/"C]. Furthermore, the ^^C^O/^^C^^^O ratio would provide 
a measure of whether the same behaviour occurs for [^''0/^**0], as 
might be expected if infall of less nuclear processed material is 
occurring. This is because the [^^0/^*0] gradient has also been 
found to incr ease with distance from the centre of the Galaxy 
( IWilsorJl999h , similarly to that of the [^^C/^^C] ratio. On the other 
hand, if the [^'^C/^^C] variation arises due to another mechanism 
(e.g. differential chemical fractionation), then the [^''O/^^O] might 
not be found to increase with increasing velocity. 



5 SUMMARY 

We have mapped a 2.5° x 0.5° region of the centre of the Galaxy 
using the Mopra radio telescope in 18 molecular lines emitting 
from 85 to 93 GHz. This incorporates most of the region known 
as the Central Molecular Zone (CMZ). The molecular maps have 
~ 40 arcsec spatial resolution and 2 km s~^ spectral resolution, 
with emission extending to Vlsr ~ ±220 km s~^. For eight 
species the emission is particularly strong and is widespread across 
the CMZ: HCN, HNC, HCO+, HNCO, N2H+, SiO, CH3CN and 
HC3N. For the other molecules mapped the emission is generally 
confined to the bright dust cores around Sgr A, Sgr B2, Sgr C 
and G1.3. This includes the isotopologues H^^CN, HN"C and 
H"C0+, as well as C-C3H2, CH3CCH, HOCO+, SO and C2H. 

As seen in several other studies of the CMZ, the molecu- 
lar emission in all these species, while widespread, is asymmetric 
about the centre of the galaxy, with roughly three-quarters com- 
ing from positive longitudes and one-quarter from negative longi- 
tudes. 35-45 percent of the total molecular emission also arises 
from the vicinity of the four brightest dust cores. The overall emis- 
sion morphology is also remarkably similar between all molecules. 
Line profiles are both very wide and complex, and do vary con- 
siderably across the CMZ. While turbulence contributes to the ~ 
30 — 50 km s^^ width of the profiles, there are also clearly several 
components along many of the sight lines. Absorption from colder, 
foreground gas is also apparent at three velocities (—52, —28 and 
—3 km s~^) for several locations and most strongly in HCN, HNC, 



and HCO - attributable to the sight line passing through spiral 
arms. 

To quantify the overall emission morphology, and its varia- 
tion between molecules, we conducted a principal component anal- 
ysis (PCA) of the integrated emission from 8 brightest species (as 
above). The first 4 components are significant, with the first compo- 
nent dominant, contributing 83 percent of the total variance mea- 
sured. It essentially provides an average integrated emission map 
of the CMZ and demonstrates that these lines are indeed similar 
in overall appearance. The second component contributes 1 1 per- 
cent to the total variance and highlights some differences in the 
bright cores of Sgr A and Sgr B2 between the HCN, HNC and 
HCO^ lines and the other species. This is attributable to these three 
lines being optically thick in these regions (which we later quan- 
tify through an optical depth analysis). The third and fourth com- 
ponents (3 and 2 percent of the variance) highlight smaller differ- 
ences between SiO, HNCO and CH3CN, reflecting the relative im- 
portance of shocks and hot molecular core emission between some 
locations. 

We have selected apertures around the bright dust cores, as 
well as for the total region mapped, in order to study line ratio 
variations and to calculate optical depths so that column densities 
and molecule masses may be determined. Dense cores species like 
CH3CN, HNCO, HC3N and N2H+ are also found to be more vari- 
able with position, compared to HCN, than are HCO^ and HNC. 

We also compare the line intensity ratios of HCN, HNC and 
HCO^ to the same ratios measured for their ^^C isotopologues. 
This has application to the interpretation of these line ratios mea- 
sured in extra-galactic sources for determining the excitation condi- 
tions, where the spatial structure remains unresolved. Of note is that 
the HCO+/HCN ratio (~ 0.6) is found to typically be 2-3 times 
the ratio of H^^C0+/H"CN. This is as a result of different opti- 
cal depth in the two lines in the ^^C ratio and indicates that with- 
out measurement of the isotopologues the diagnostic ratio obtained 
may be in error. However, in all cases for the CMZ the diagnostic 
ratios are found to lie well away from the XDR-excitation region, 
and within the PDR region - as expected given our knowledge of 
the overall molecular environment in the centre of the Galaxy. 

An optical depth analysis using the '^^C to '^^C isotopologue 
line ratios finds modest optical depth for HNC and HCN, with 
r ~ 1 for HCO^. Correcting the optical depth allows column 
densities and molecular masses to be determined for the CMZ. 
Around 50 Mq of [HCN+HNC+HCO+] exists with the region, for 
instance. 

The ratios of the ^^C to ^"^C isotopologues are found to in- 
crease with red-shifted velocity in several locations, indicating a 
decreasing optical depth in the line wing. This is particularly no- 
ticeable for HCO^. However at the extreme red-end {Vlsr > 
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+ 100 km s^^) the ratio exceeds 24, the value determined for the 
12q^13q isotope r atio in the centre of the Galaxy. This result has 
also been found bv lRiguelme et al. I( l2010bh for other selected loca- 
tions in the centre of the Galaxy. We discuss whether this provides 
evidence for infall of less nuclear-processed gas (where the isotope 
ratio is higher) onto the xi and X2 orbits in the central regions, so 
refuehng the gas reservoir in the CMZ. 

Line luminosities, relative to that of CO, vary from ~ 10 per- 
cent for HCN to ~ 1 percent for SiO. The luminosities are also typ- 
ically 0.1 — 10 percent of the corresponding values that have been 
measured in other galaxies, though this comparison is biased by the 
selection of the external galaxies studied, as well as the comparison 
of the CMZ to regions several kpc in extent in those galaxies. How- 
ever, the relative line luminosities (to CO) are not dissimilar, with 
the CMZ lying towards the lower end of the range in HCN/CO ra- 
tios observed. This might be attributable to a somewhat lower star 
formation rate, related to the availability of dense gas, than in the 
comparison galaxies. 

The full data set, comprising the data cubes for the 20 emis- 
sion lines, is publicly available for further analysis and may be 
obtained by contacting the authors. It will also be made avail- 
able through the Australia Telescope National Facihty archives 
(www.atnf.csiro.au). A further data set, comprising 24 lines mea- 
sured in the 7 mm band between 42-50 GHz with the Mopra tele- 
scope, is also being prepared for publication. It includes several 
lines from lower energy levels of the same molecules as presented 
here, so allowing an excitation analysis to be conducted. 
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